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This contribution describes the HOPPET Fortran 95 package for carrying out DGLAP evo- 
lution and other common manipulations of PDFs. The PDFs are represented on a grid in 
x-space so as to avoid limitations on the functional form of input distributions. Good speed 
and accuracy are obtained through the representation of splitting functions in terms of their 
00 ' convolution with a set of piecewise polynomial basis functions. We briefly describe the struc- 

f^ , ture of the code and discuss its quantitative performance. Finally we comment of future 

^3 ' directions for the program's development. 

(N , 

Introduction There has been considerable discussion over the past years of numerical solutions 
of the DGLAP equation for the QCD evohition of parton distributions (PDFs), 
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'11^. There exist two main classes of approaches: those that solve the equation directly in a;-space and 

^L{ those that solve it for Mellin transforms of the parton densities, defined as 
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^ ' and subsequently invert the transform back to x-space. iV— space based methods are of interest 

OO , because the Mellin transform converts the convolution of eq. ^ into a multiplication, resulting 

in a continuum of independent matrix differential (rather than integro-differential) equations, one 
for each value of N, making the evolution more efficient numerically. 

The drawback of the Mellin method is that one needs to know the Mellin transforms of both the 
splitting functions and the initial conditions. The x-space method is in contrast more flexible, since 
the inputs are only required in x-space; however it is generally considered less efficient numerically, 
(^ , because of the need to carry out the convolution in eq. H]). 

Despite it being more difficult to obtain high accuracy with x-space methods, their greater 
flexibility means that they are widespread, serving as the basis of programs like QCDNUM [2] 
^ ■ or CANDIA [3], and used also by the CTEQ and MSTW global fitting collaborations. A novel 

approach which combines advantages of the iV— space and a;— space methods is discussed in [4]. 

HOPPET [5], the program described in this contribution [1], uses higher-order methods both 
for the x-integrations and Q evolution. Apart from its mainstream PDF evolution capabilities, it 
also provides access to a range of low and medium-level operations on PDFs which should allow 
a user to extend the facilities already provided. 

The functionality of hoppet has been present in its predecessors for several years, and a first 
public release with complete documentation was recently presented [5]. Those predecessors have 
been used in a number of different contexts, like resummation of event shapes in DIS [6], automated 
resummation of event shapes [7] , studies of resummation in the small- a; limit and in a posteriori 
inclusion of PDFs in NLO final-state calculations [8] , as well as used for benchmark comparisons 
with Pegasus [5] in [TU] . 

Program structure The numerical techniques upon which HOPPET is based are presented in [5]. 
Let us only recall here that it uses higher-order methods both for the x-integrations and Q evo- 
lution. It combines this with multiple grids in x-space: a high-density grid at large x where it 
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STREAMLINED INTERFACE 



METHOD 


DESCRIPTION 


Initialisation 




hoppetStart (dy ,nloop) 


Sots up a compound grid with spacing in In Xjx of dy at small a:, 

extending to y — 12 and numerical order = — 5. 

The Q range for the tabulation will be 1 GeV < Q < 28 TcV, 

dlnlnQ=dy/4 and the factorisation scheme Is MS 


hoppetStartExtended (ymax , dy , Qmin , 
Qmax , dlnlnQ ,nloop , order , fact scheme ) 


More general initialisation 


hoppetSetFFN(f ixed_nf ) 
hoppetSetVFN (mc , mb, mt) 


Set heavy flavour scheme 


alphas = hoppet Alphas (Q) 


Accessing the coupling 


Normal evolution 




hoppetEvolve (asQ , QOalphas , 
nloop , muR_Q , LHAsub , QOpdf ) 


PDF evolution: specifics the coupling asQ at a scale QOalphas, 

the number of loops for evol., nloop, 

the ratio (muR_C]) of ren. to fact, scales. 

the name of a subroutine LHAsub with an LHAPDF-llkc Interface 

and the scale QOpdf at which one starts the PDF evolution 

Note: LHAsub only called at scale QOpdf 


hoppetEval(x,Cl,f ) 


On return, f (-6:6) contains all flavours of the PDF set 
(multiplied by x) at the given x and Q values 


Cached evolution 




hoppetPreEvolve (asQ , QOalphas , 
nloop, muR.Q, QOpdf) 


Preparation of the cached evolution 


hoppetCachedEvolve (LHAsub) 


Perform cached evolution with the initial condition 

at QOpdf from a routine LHAsub with LHAPDF-llkc Interface 

Notice LHAsub only called at scale QOpdf 


hoppetEval(x,Q,f ) 


On return, f (-6:6) contains all flavours of the PDF set 
(multiplied by a:) at the given x and Q values 
[as for normal evolution] 



Table 1: Reference guide for the streamlined interface. 



is hardest to obtain good accuracy, and coarser grids at smaller x where the smoothness of the 
PDFs facilitates the integrations. One of the other crucial features of the program is that it pre- 
calculates as much information as possible, so as to reduce the evolution of a new PDF initial 
condition to a modest set of addition and multiplication operations. 

The functionalities in hoppet can be accessed from two interfaces. The first one is a general 
interface, available in F95 only, required for those who want to delve into the innards of the 
program, and which provides access to a range of low and medium-level operations allow to 
extend hoppet basic capabilities. The essential hoppet functionalities are also accessible from a 
streamlined interface, which can be accessed from F97, F95 and C++. 

A summary of the most important modules in the streamlined interface can be seen in Table 
[TJ It is clear that, thanks to the structure of this interface, the evolution of PDFs with a large 
variety of options can be performed with few selected routines. The analog reference guide for the 
general interface can be found in [S] . 

Often one simply wishes to provide PDFs at some initial scale and then subsequently be able 
to access it at arbitrary values of x and Q. For this purpose it is useful (and most efficient) in 
HOPPET to produce a table of the PDFs as a function of Q^, which then allows for access to the 
PDFs at arbitrary x, Q using an interpolation. 

It is also possible to prepare an evolution in cached form. This is useful if one needs to evolve 
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Tabic 2: Contributions to the run time in eqs. ([3]) for dy = 0.2 and dlnlnQ = 0.05 and standard 
values for the other parameters (on a 3.4GHz Pentium IV (D) with 2 MB cache). 



many different PDF sets with the same evolution properties (coupling, initial scale, etc.), as is the 
usual situation in global analyses of PDFs, the initialisation may take a bit longer than a normal 
evolution (2-10 times depending on the perturbative order), however, once it is done, cached 
evolutions run 3-4 faster than a normal evolution. 



Performance hoppet's correctness has been established with a reasonable degree of confidence 
in the benchmark tests [TD] where it was compared with the Mellin space based evolution code 
QCD-Pegasus H]. 

For use in most physical applications, an accuracy in the range 10^"^ to lO^'* is generally more 
than adequate. The critical issue in such cases is more likely to be the speed of the code, for 
example in PDF fitting applications, moppet's accuracy and speed both depend on the choice of 
grid (in y = Inl/a;) and the evolution and/or tabulation steps in Q. 

Fig. [T] shows the relative accuracy e as a function of x for dy = 0.05 and dlnlnQ = dy/4. 
The relative accuracy e is poorest as one approaches x = 1, where the PDFs all go to zero very 
rapidly and so have divergent logarithmic derivatives in x, adversely affecting the accuracy of the 
convolutions. This region is always the most difficult in x-space methods, however the use of 
multiple subgrids in x allows to one to obtain good results for x < 0.9 which is likely to be the 
largest value of any phenomenological relevance. 

The time spent in moppet for a given analysis can expressed as follows, according to whether 
or not one carries out cached pre-evolution: 



'^no prc-cv ^s ^" ^a^a i rii\Zi -f- ^xQ ^xQ ) 7 
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where ts is the time for setting up the splitting functions, ria is the number of different running 
couplings that one has, ta is the time for initialising the coupling, rii is the number of PDF initial 
conditions that one wishes to consider, ti is the time to carry out the tabulation for a single initial 
condition, Hxq is the number of points in x, Q at which one evaluates the full set of flavours once 
per PDF initial condition; in the case with pre-prepared cached evolution, tp is the time for a 
preparing a cached evolution and tc is the time for performing the cached evolution. Finally txQ 
is the time it takes to evaluate the PDFs at a given value of {x, Q^) once the tabulation has been 
performed. 

The various contributions to the run-time are shown in Table[2]for dy — 0.2 and dlnlnQ = 0.05 
(giving an accuracy ~ 10^"*), for various compilers. In a typical analysis where run-times matter, 
such as a PDF fit, it is to be expected that the time will be dominated by tc (or ti). However, in 
the typical case of global PDF fits, for which the number of x, Q points is rather large (> 3000), 
it will be UxQtxQ that takes the most timco. 

We study tc and ti in more detail in Fig. [l] where we relate them to the accuracy obtained 
from the evolution. As one would expect, studying just the 'guds' flavours for x < 0.7 one obtains 



^ With these numbers, it is easy to check that a global fit with nj;Q 
in less than half an hour. 



3000 and Ui ~ 10'' could be completed 
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Figure 1: Left: the relative accuracy e of the least well determined flavour channel at each x, Q 
point, shown as a function of x for many Q values. Right: the same quantity obtained as a function 
of the time taken to perform the evolution. 

better accuracy for a given speed than with all flavours for x < 0.9. Overall one can obtain 10^^ 
accuracy with tc — 10""^ s and 10~^ accuracy with tc — 10^^ s. 

The timings shown here are roughly similar, for accuracies ~ 10^^, to those obtained with 
the A^-space code Pegasus [5] when the number of x, Q points to be evaluated is O (10'^). To the 
best of our knowledge, other NNLO evolution codes published in recent years seem generally less 
competitive either in terms of accuracy or speed. 

Future directions HOPPET is an x-space evolution code that is novel both in terms of the accuracy 
and speed that it provides compared to other x-spacc codes, and in terms of its interface, designed 
to provide a straightforward and physical way of manipulating PDFs beyond the built-in task of 
DGLAP evolution. 

Work in progress for future releases include the computation of deep-inelastic scattering struc- 
ture functions and reduced cross sections. Furthermore, a physical feature absent from mainstream 
PDF evolution codes is that of evolution that includes matching with various types of resummed 
calculations. There is work in progress to implement the ABF and CCSS small- a; resummations 
[IT] into HOPPET through the ability to manipulate general, user-defined, interpolated splitting 
functions and coefficient functions. 
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